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Effect of radiolysis of TODGA on the extraction of TODGA/n- 


dodecane toward Eu(III): An experimental and DFT study 


Abstract: N,N,N’,N’-Tetraoctyl diglycolamide (TODGA) is one of the most 
promising extractants tailored for high-level liquid radioactive waste treatment 
during nuclear fuel reprocessing. The y-radiolysis of TODGA (0.2 mol/L) in n- 
dodecane (nDD) solution with and without pre-equilibrated 3.0 mol/L HNO; was 
investigated using HPLC and UPLC-QTOF-MS and compared with the y- 
radiolysis of neat TODGA in this study. With increased absorbed doses, the 
concentration of TODGA decreased exponentially for the studied systems. 
Moreover, pre-equilibration with HNO; (3.0 mol/L) slightly influenced the y- 
radiolysis of TODGA in nDD. Seven radiolytic products generated from the 
rupture of the C — C, C — O, and C —N bonds in TODGA were identified in the 
studied extraction system. The influence of y-radiation on TODGA/nDD for the 
extraction of Eu(IIJ) was evaluated using the first combination of extraction 
experiments and density functional theory (DFT) calculations, in which the 
complexations of Eu(IIJ) with TODGA and its radiolytic products were 
systematically compared. Based on the radiolysis kinetic model of TODGA, the 
slope curve of the distribution ratio of Eu(III) (Dz,,) and the absorbed dose, and 
fluorescence titration analysis, the empirical equation of the absorbed dose and 
Dry, Was obtained successfully. Below 300 kGy, the experimental Dgu agreed 
well with the obtained empirical equation for TODGA/nDD. Conversely, at a 
high absorbed dose, the experimental Dg, was higher than the theoretical Dz, 
based on the empirical equation because the radiolytic products of TODGA with 
similar coordination structures still possessed partial complexation toward Eu(II), 
which was confirmed by DFT calculations. This work provides a method to 
predict the extraction distribution ratio of an irradiated extractant system and to 


understand the complex extraction process. 
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1 Introduction 


Spent nuclear fuel (SNF) reprocessing generates high-level liquid radioactive 
waste (HLLW), which is composed of unextracted U and Pu, long-lived fission 
products (11I, Tc, !°°Cs, and Zr), short-lived fission products (?°Sr and !?7Cs), minor 
actinides (MAs) (Np, Am, and Cm), and stable lanthanide isotopes (Eu, Nd, La, Tb, Pr, 
Gd, Sm, Ce, and Pm)[1—10]. To reduce the radiotoxicity and volume of HLLW and 
improve resource utilization, partitioning and transmutation strategies have been 
explored to shorten the half-life of long-lived radioisotopes. However, the transmutation 
efficiency of MAs is limited by the high thermal neutron capture cross-sections of 
several lanthanide isotopes owing to competition for available neutrons[11, 12]. Thus, 
separating MAs from lanthanides (Lns) before transmutation is crucial. 

During the liquid—liquid extraction process, N,N,N’,N’-tetraoctyl diglycolamide 
(TODGA) has proven to be efficient in extracting actinides (Ans) and lanthanides [2, 
13-16]. Furthermore, TODGA is a tridentate ligand that shows high distribution ratios 
for Lns(IID)/Ans(III) because the hydrophilic groups exhibit stable metal coordination in 
the aqueous phase and the hydrophobic groups ensure its solubility in the organic 
phase[17, 18]. It has been used in some processes, such as the innovative Selective 
ActiNide EXtraction (i-SANEX) process[19, 20] and the EUROpean Grouped Actinide 
EXtraction (EURO-GANEX) process[21, 22]. The i-SANEX process was first arranged 
utilizing an organic solution with TODGA (0.2 mol/L) and 1-octanol (5% v/v) to enable 
co-extraction of trivalent Lns and Ans from simulated HLLW. At the same time, Pd and 
Zr were masked by cyclohexane diamine tetraacetic acid (CDTA). Second, the co- 
extracted Mo and Sr were stripped using 3.0 mol/L HNO3 aqueous solution with CDTA 
(0.05 mol/L) and oxalic acid (0.2 mol/L) as scrubbing solvents, followed by a back- 


extraction step using an aqueous solution containing 2,6-bis(5,6-di(sulfophenyl)-1,2,4- 


triazin-3-yl)pyridine (SO3-Ph-BTP) (1.8107 mol/L) in HNO; (0.35 mol/L) for 
selective stripping of Ans(III). The final step was the stripping of Lns(IID) using a citrate 
buffer solution. The experimental results obtained by Wilden et al.[19] demonstrated the 
effectiveness of this process. In the Ans product fraction, Cm(III) and Am(III) were 
efficiently stripped (more than 99.8%) and regained within six stages. The purity of the 
Cm(III) and Am(III) products was found to be very good, with only Ru (0.4%), Sr 
(0.3%), and Lns(III) (0.1%) contamination. Lns(III) stripping is efficient (more than 
99.5%) within only four stages using a stripping solution containing citric acid with a 
pH of 3. The EURO-GANEX process consists of two cycles. In the first cycle, 
hydrogenated tetrapropene (TPH) containing N,N-di (2-ethylhexyl) isobutyramide 
(DEHiBA) (1.0 mol/L) was applied to extract U(VI) selectively. In the second cycle, the 
extraction of Cm(II), AmdII), Pu(dlV), Np(VD, and Lns(III) from the first cycle 
raffinate was achieved using N,N’-dimethyl-N,N’-dioctyl-2(2-hexyloxyethyl)- 
malonamide (DMDOHEMA) (0.5 mol/L) and TODGA (0.2 mol/L) in a kerosene 
diluent, while CDTA (0.055 mol/L) was applied as a masking agent, to prevent the 
extraction of Zr(IV) and Pd(II). The Sr and Fe were extracted using HNO; (0.5 mol/L) 
from the loaded solvent, and trans-uranium elements (TRU) were back-extracted by 
applying an aqueous solution containing aceto-hydroxamic acid (AHA) and SO3-Ph- 
BTP, while Lns(III) was retained in the organic solution. The EURO-GANEX process 
was successfully demonstrated by the hot test performed by Malmbeck et al.[23]. 
Moreover, Pu([V), Np(VI), AmdII), and Lns(III) were efficiently co-extracted from the 
raffinate of the EURO-GANEX process during the first cycle; Mo and Zr were also co- 
extracted, whereas the other fission products were successfully refused. Within six strip 


stages, 99.91% Np, 99.85% Pu, 99.90% Am, and co-extracted Mo and Zr were well 


back-extracted, whereas 99.9% Lns(III) was retained in the organic solution within four 
strip stages using an aqueous solution containing AHA and SO3-Ph-BTP. 

Previous studies have shown that TODGA provides significantly higher 
distribution ratios in alkanes such as n-dodecane (nDD) than in other studied diluents, 
such as secondary ketones, primary alcohols, and secondary alcohols[24, 25]. The 
extraction system of TODGA/nDD can be applied to separate MAs and Lns from other 
fission products. Once separated, the extraction of individual elements can be simplified 
with additional processing, culminating in the separation of MAs from Lns [26-28]. 

Because the extraction system is exposed to the radioactive environment formed by 
the radionuclides, its radiation stability is an important factor to be examined before its 
practical application[29]. In the solvent extraction process, the y-radiolysis of TODGA 
is problematic because of the highly radioactive HLLW. Because the loss of TODGA 
and the radiolytic products of TODGA formed during y-radiation can adversely affect 
the extraction performance in the liquid—liquid extraction process toward metal ions, it 
is essential to measure the y-radiolysis yield and radiolytic products of TODGA and to 
study their effects on the performance of extracting metal ions. Radiolysis of TODGA 
in nDD under different conditions has been investigated in some studies [30, 31]. It was 
found that the magnitude of a-radiolysis of TODGA is lower than that of y-radiolysis, 
and they have similar radiolytic products[32]. This is attributed to the fact that the linear 
energy transference of a particles is higher than that of y particles, and the radical yield 
of nDD decreases by recombination in a-radiolysis. Although the y-radiolysis and 
extraction of TODGA have been investigated, the effect of y-radiolysis and radiolytic 
TOGDA products on extraction in the TODGA/nDD system has rarely been reported. In 
addition, if we could predict the extraction performance of the irradiated extractant 


system using equations, it would be very beneficial and important in controlling the 


partition process. To the best of our knowledge, the mathematical correlation between 
the radiolysis of the extractant and the distribution ratio has seldom been studied. 

In this work, because a and y rays result in similar radiolysis products of TODGA 
and y-radiation is easier to carry out in the laboratory, the radiolytic behavior and y- 
radiation stability of TODGA in the 0.2 mol/L TODGA/nDD system (TODGA/nDD) 
and the 0.2 mol/L TODGA/nDD that was pre-equilibrated with 3.0 mol/L HNO; system 
(HNO3-TODGA/nDD) were studied qualitatively and quantitatively. They were also 
compared with radiolysis of neat TODGA. The content and radiolysis products of 
TODGA were determined using HPLC and UPLC-QTOF-MS, respectively. 
Considering the chemical similarity between trivalent lanthanides and minor actinides, 
and to avoid the dangers of handling radioactive isotopes in standard laboratories, Eul) 
was selected as a model metal ion in our extraction experiment to assess the extraction 
performance of TODGA/nDD after y-radiation. The density functional theory (DFT) 
calculation, which can estimate the coordination abilities of TODGA and its radiolytic 
products, was applied to explain the effect of the radiolytic products on our extraction of 
Eu(iI]). For the first time, a combination of theoretical calculations and extraction 
experiments was used to demonstrate that the radiolytic products maintain partial 
complexation, which capably explains the extraction behavior of irradiated samples. For 
the TODGA/nDD system, an in-depth understanding of the effects of y-irradiation on 


the extraction behavior was obtained from this study. 


2 Experimental 


Materials 


Sichuan University synthesized TODGA (298%); its purity was confirmed by FT- 
MS (ESI*) and 'H NMR (Figure S1 — 2). Moreover, Eu(NO3)3-6H20 (99.99%) was 
obtained from MREDA Technology Co. Ltd. Formic acid, liquid chromatography-mass 
spectrometry (LC-MS) grade, was purchased from Thermo Fisher Scientific. Methanol 
(LC-MS grade) was purchased from Honeywell Trading (Shanghai) Co., Ltd. The 
aqueous dilution was performed using ultrapure water (18.2 MQ-cm). All other 
chemicals were of analytical grade, and the chemicals were not further purified. 


Gamma irradiation 


TODGA and TODGA/nDD solutions equilibrated with and without HNO; (3.0 
mol/L) before y-radiation were irradiated in the air at room temperature (25 + 5 °C) 
using a £Co source (Department of Applied Chemistry in the College of Chemistry, 
Peking University, China) with an absorbed dose rate of 6.6 kGy/h, as determined by a 
Fricke dosimeter. 


Quantitative analysis of TODGA after irradiation 


An HPLC method with UV detection (SPD-16, Shimadzu, Japan) was used to 
quantitatively analyze the TODGA concentration in the irradiated samples. Additionally, 
HPLC was carried out at 25 + 5 °C using an Agilent HC-C18 liquid chromatography 
column (5 um, 4.6 mm x 150 mm) to achieve chromatographic separation with a UV 
wavelength fixed at 210 nm. The aqueous component was ultrapure water with formic 
acid (1%o v/v) (A); the organic component was methanol (B). The elution program in 
our liquid chromatography was listed hereafter: 0-10 min, 30-100% B; 10-20 min, 100% 
B; 20-22 min, 100-30% B. The error values in the quantitative analysis experiments 


were within 5%. The flow rate was 1.0 mL/min, and the injection volume was 20 uL. 


Identification and semi-quantification of radiolytic products 


Ultra-high pressure liquid chromatography (UPLC) (ACQUITY I-Class, Waters, 
USA) was carried out using an Agilent HC-C18 liquid chromatography column (5 um, 
4.6 mm x 150 mm) at 40 °C to achieve chromatographic separation. The aqueous 
component was ultrapure water with formic acid (1%o0 v/v) (A), and the organic 
component was methanol with formic acid (1%o v/v) (B). The liquid chromatography 
elution program was listed hereafter: 0—4 min: 30-100% B; 4-13 min: 100% B; 13- 
13.1 min: 100-30% B; 13.1—15 min: 30% B. The flow rate was 0.8 mL/min with a 1 uL 
injection volume. The mass-spectrometer conditions of the quadrupole time-of-flight 
mass spectrometry (QTOF-MS) (Vion IMS QToF, Waters, USA) are listed here: 
desolvation temperature: 300 °C; source temperature: 120 °C; desolvation gas flow: 800 
L/h; capillary voltage: 2.8 kV; cone gas flow: 50 L/h; positive mode. The data were 
acquired and processed using a Waters UNIFI Scientific Information System with a 
mass target match tolerance of less than 5 ppm. 


Extraction experiment 


The TODGA was dissolved with nDD to obtain a 0.7 ml organic phase with a 
concentration of 0.2 mol/L. The Eu(NO3)3*6H20 was dissolved in 3.0 mol/L HNO; to 
obtain a 0.7 ml aqueous phase with a concentration of 1000 mg/L. Pre-equilibration of 
the organic extraction phase with 3.0 mol/L HNO; was performed before extraction. A 
vortex mixer at a speed of 2500 rpm (LPD2500, Leopard Scientific Instrument Co., Ltd, 
Beijing) was used to perform the extraction experiments at 298 + 1 K. To achieve 
complete separation between the two phases, the mixtures were centrifuged for 1 min at 
5000 rpm in a centrifuge (TGL-16M, Cence, Hunan). After dilution with ultrapure 
water, the aqueous phase was detected using an ICP-MS (ELEMENT XR, Thermo 


Scientific, USA) to estimate the concentration of Eu(II). 


The distribution ratio of Eu(I) (Dgu) was calculated using Eq. (1). The phase ratio 
equals 1:1 in the extraction experiment. 

Dey = (Ci — Cf) /C (1) 
where Cy and C; represent the final and initial Eu (III) concentrations in the diluted 
aqueous phase, respectively. The extraction had an error rate of less than 5%. 
Fluorescence spectra titration 


The fluorescence emission spectra of the titration experiments were investigated at 
298 + 1 K in a cuvette cell with a 1 cm path length using a fluorescence 
spectrophotometer (F7000, Hitachi, Japan). The concentration of TODGA was 0.1 
mol/L in acetonitrile and the initial concentration of Eu(II]) was 0.001 mol/L in 2.0 mL 
acetonitrile. Next, 4 uL of TODGA solution was added at a time. The solution was 
mixed using a vortexer (VORTEX 3, IKA, Germany) for 5 min after the addition of 
each ligand. The titration experiments used an excitation wavelength of 395 nm; the 
bandwidth of the excitation and emission was 2.5 nm. The interval was 1 nm and the 
spectra were recorded between 550 and 700 nm. Single-component spectra of the 
metal—ligand complexes and metal solvent species were obtained using the HyperSpec 
program. 


Theoretical calculations 


DFT calculations incorporating electron correlation effects were performed at the 
level of B3LYP using the Gaussian 09 package [33-35]. The relativistic effects of Eu 
atoms were investigated using the quasi-relativistic effective core potentials (RECPs) 
and associated valence basis sets, which were improved by the Stuttgart and Dresden 
groups[36—40]. The large-core RECPs used in the structural optimization of Eu(IID) 
contained 52 electrons[39, 40]. All other C, N, O, and H atoms utilized the 6-31G (d) 


basis set. Geometrical optimization and the electronic calculations of all structures were 


first performed at the level of B3LYP/6-31G(d)/RECP under the gas phase conditions. 
The Gibbs free energy (Gg), entropy (Sg), and enthalpy (Hg) were obtained under the gas 
phase conditions at 298.15 K using the same theory level. To obtain the Gibbs free 
energy (Gaq), entropy (Saq), and enthalpy (Haq) of the species in the nDD phase at 
298.15 K, all species were optimized at the B3LYP/6-31 G (d)/RECP level in nDD to 
better predict their solvation energy[41]. This theory is based on the SMD universal 


continuum solvation model[42]. 


3 Results and discussion 


Quantitative analysis of TODGA under different conditions after y-ray 
irradiation 


TODGA/nDD, HNO3-TODGA/nDD, and neat TODGA were subjected to y- 
irradiation with absorbed doses ranging from 20 to 500 kGy. The irradiated samples 
were measured quantitatively using HPLC-UV and their radiolytic stability was 
analyzed. Figure 1 shows that the TODGA concentration decreased exponentially with 
increasing absorbed doses, suggesting pseudo-first-order degradation kinetics for 
TODGA/nDD, HNO3-TODGA/nDD, and neat TODGA. By comparing TODGA/nDD 
and HNO3-TODGA/nDD, it was found that HNOs had little influence on the radiolysis 
of TODGA. Moreover, in the nDD solution, the radiolysis rate of TODGA was higher 
than that of the neat TODGA. A dodecane-induced "sensitization effect" is responsible 
for the sensitization of TODGA to radiolysis[43] due to the fact that the radiolysis of 
TODGA is accelerated by nDD solvents because of the efficiency of nDD in 


transferring positive charges [44-46]. 


Because the radiolysis behavior of TODGA depends on its concentration in 
pseudo-first-order kinetics, the rate of radiolysis causes an exponential drop in the 


TODGA concentration. The radiolysis rate equation is defined by Equation (2): 


dc 
= A = Cropga ' K (2) 


The solution to the rate equation is presented as Equation (3): 


CTODGA _ „-—kR 
Tone = e 3) 
TODGA 


where Cropga (mol/L) is the concentration of TODGA, Cfopga (mol/L) is the 
original concentration of TODGA before y-radiation, R (kGy) is the absorbed dose, and 
the dose constant k (kGy”') is the coefficient obtained by fitting the exponential curve. 

The dose constant k can be used to scale the absorbed dose when exponential 
behavior occurs[47]. Using all the information acquired during the irradiation 
experiment provides more accurate results, and standard statistical methods can be used 
to assess the accuracy of the results. The application of k to characterize an 
exponential curve in terms of the absorbed dose is very similar to that of traditional 
time-dependent kinetics. As expressed by exponential Equation (3), the reciprocal 


absorbed dose is governed by the first-order rate law instead of the reciprocal time. 


e Neat TODGA 
4 TODGA/nDD 
m HNO,-TODGA/nDD 


y = exp(-6.54 x 10x) 
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Figure 1 Crgpga/Cropca of irradiated TODGA/nDD, HNO3-TODGA/nDD, and neat 
TODGA as a function of the absorbed dose (R) (y: Cropga/Cfopga »X: R) 


If the dose constants k of different systems can be obtained, Equation (4) can be 

used to calculate the G value (umol/J), which is known as the radiation chemical yield: 

G = 10ëCopgakp ~! (4) 

where p is the solution density (kg/L) (Prapga= 9.91 and Po.2 mol/L ropGa/npp= 9.77 
(25 °C). The measurement results are shown in Table S1. 

Table 1 shows the dose constants k and radiation chemical yields G of irradiated 
TODGA/nDD, HNO3-TODGA/nDD, and neat TODGA. The dose constant k of 
TODGA in nDD is higher than that of neat TODGA because of the “sensitization effect,” 
while the radiation chemical yields G of TODGA/nDD and HNO3- TODGA/nDD are 


lower than that of neat TODGA because of the different initial concentrations of 


irradiated TODGA. These results agree with the literature in which Sugo et al. reported 
that the k value of 0.2 mol/L TODGA/nDD was 2.2x10° kGy"'[43, 48]. By using the 
k value, we can predict that the absorbed doses for the half-loss of TODGA 
concentration (Ros) in irradiated TODGA/nDD, HNO3-TODGA/nDD, and neat 


TODGA will be approximately 330, 315, and 1060 kGy, respectively. 


Table 1 Dose constants k, radiation chemical yields G, and absorbed dose for half-loss 


of concentration of TODGA Ro.s of irradiated TODGA/nDD, HNO3-TODGA/nDD, and 


neat TODGA 
Samples CRopga mol/L) k (kGy"') G (umol/J)  Ro.s (kGy) 
TODGA/nDD 0.20 (2.08+0.05)x10° 0.540.013 330 
HNO3-TODGA/nDD 0.20 (2.20+0.03)x10° 0.570.008 315 
Neat TODGA 1.57 (6.54+0.31)x107 1.130.053 1060 


Identification and semi-quantitative analysis of radiolytic products 


Ultra-high-performance liquid chromatography with quadrupole time-of-flight 
mass spectrometry (UPLC-QTOF-MS) can be used to identify and quantify complex 
mixtures of unknown compounds owing to its high resolution and mass accuracy. 

The major radiolytic products of TODGA were identified using QTOF-MS; these 
are presented in Table 2, and their UPLC-QTOF-MS spectra are shown in Figures S39. 
Seven mass signals were identified and attributed to the possible radiolytic products of 
TODGA (P1— P7, Table 2). These radiolytic products were formed by the rupture of the 
C—N, C—O, and C —C bonds of TODGA. Four radiolysis routes after irradiation are 


shown in Scheme 1. 


Table 2 Radiolytic products of TODGA determined by the QTOF-MS 
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Scheme 1 Radiolysis routes of TODGA after y-irradiation 


Figure 2 shows the semi-quantitative analysis of the radiolytic products Pı— P7 in 
the irradiated HNO3-TODGA/nDD and neat TODGA. The absorbed doses ranged from 
20 to 500 kGy. Seven radiolytic products were found to exist in the two systems. The 
detector counts of the radiolytic products Pı — P7 roughly increased with increasing 
absorbed dose. However, some radiolytic products, such as Pi in irradiated neat 
TODGA, showed a decrease in detector counts with absorbed doses from 400 kGy to 
500 kGy. This is probably because the contents of the radiolytic products during 


irradiation increased from the radiolysis of the TODGA and decreased by self-radiolysis. 
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The concentration’s rate of increase was less than the rate of decrease with absorbed 
doses from 400 to 500 kGy for Pi in irradiated neat TODGA. The detector counts of all 
radiolytic products in irradiated TODGA/nDD were higher than those in neat TODGA 
(Figure S10), indicating the “sensitization effect” of nDD. This agreed well with the 
above result of the radiolysis kinetics of the TODGA with or without nDD. However, 
the relative detector counts of some radiolytic products differed between the two 
irradiation systems. For example, the detector count of Po was less than that of P4 in 
irradiated neat TODGA; however, the converse was true for irradiated HNO3- 
TODGA/nDD. This can be attributed to the different reaction mechanisms, wherein 
irradiated neat TODGA produces radical cations by direct ionization reaction. 
Conversely, in irradiated TODGA/nDD, nDD radical cations may transfer their charge 


to TODGA molecules, leading to radiolysis of TODGA[43, 49]. 
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Figure 2 Semi-quantitative analysis of radiolytic products P1 — P7 in: (a) irradiated neat 
TODGA, and (b) HNO3-TODGA/nDD with absorbed doses ranging from 20 to 500 
kGy 


Influence of y-ray irradiation on the extraction of TODGA system 


toward Eu(iII) 


HNO3-TODGA/nDD and neat TODGA were subjected to y-ray irradiation at 
absorbed doses ranging from 50 to 500 kGy. The irradiated neat TODGA was diluted 
with nDD to 0.2 mol/L. These organic phases were used to extract Eu(II) from the 3.0 
mol/L HNO; aqueous phase to evaluate the extraction performance. 

As shown in Figure 3, Dz, decreased as the absorbed dose increased, but it was 
slower for irradiated TODGA than for irradiated HNO3-TODGA/nDD, indicating that 
y-ray irradiation of the ligand significantly affected the extraction of Eu(III). This 
further proves that the presence of nDD increased the radiolysis of TODGA. However, 
the HNO:-TODGA/nDD system after y-ray irradiation still possesses a high Dry 


(2.2x10°), even at 500 kGy. 
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Figure 3 Dp„ of irradiated neat TODGA and HNO3-TODGA/nDD for extraction of 
Eu(IID versus the absorbed dose. 


The loss of extractants was the main cause of the decline in extractability after 
irradiation. Owing to the presence of TODGA and its liquid radiolysis products in the 
irradiated organic extraction phase, the complexation of radiolytic products on Eu(II) 
was studied using DFT calculations. This determined the effect if any of the radiolytic 
products of TODGA on the extraction behavior of the irradiated extraction system. 


Slope analysis and fluorescence titration 


Before DFT calculations were conducted to investigate the complexation of 
TODGA and its radiolytic products with Eu(III), we assessed the stoichiometry of Eu 
ions and TODGA ligands in the extraction reaction, which could be obtained by the 
slope analysis in the curve of log Dm — log[L](org). 

The extraction reaction of Eu(II) by TODGA is expressed by Equation (5): 


(5) 


Eu(aq) + 3NO3 „4) + NTODGA org) = Eu(NO3)3TODGA 


"(org.) 


The extraction equilibrium concentration constant, Key, is described as Equation 


(6): 
[Eu(NO3)3TODGAn|(org,) 
Kex = 6 
eX [Eu3+](aq) [N03 Nag [TODGA org) (6) 
The distribution ratio, Dgu, is defined as Equation (7): 
__ [Eu(NO3)3TODGAn| (org) 


Equation (7) was substituted into Equation (6), and Equation (6) was then 
transformed into a logarithmic form to obtain Equation (8). 

log Dey, = log Kex + nlog[TODGA] (o,g) + 3 log[NO3 |(aq) (8) 

The n value can be calculated from Equation (8), which represents the average 

number of TODGA coordinated to one metal ion. The plot in Figure 4 of log Deu — 


log[TODGA] org) at a constant aqueous solution of 3.0 mol/L HNO; shows that the 


slope is 3.10 + 0.11, which indicates that Eu(III) would form complexes with three 


TODGA molecules. 
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Figure 4 Plot of log Dg, versus log[TODGA] (org) at a constant aqueous 3.0 mol/L 


HNO; solution (x: log[TODGA] org), y: log Deu) 


Complexation studies were also performed with TODGA and Eu(II) using 
fluorescence titration. Figure 5 displays the normalized fluorescence emission spectra 
with different ratios of Eu(III) and TODGA resulting from the transitions of 5Do—7F1 
and 5Do—7F2. The transition of 5Do—7F: at 593 nm, owing to the magnetic-dipole 
transition, was considered independent of the ligand field. The transition of 5Do—7F2 at 
617 nm owing to the electric—dipole transition is considered to be hypersensitive to the 
ligand field, and its intensity depends on the coordination symmetry around the metal 


ions[50]. With the addition of the TODGA ligand, the change in the ligand field led to a 
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gradual transformation of the single peak at 617 nm into a double peak, which was 
finally formed at M: L = 1:3. This indicated that the inner coordination sphere of EuD) 
was gradually occupied by TODGA with the addition of the ligand, which is similar to 
some results from the literature[51, 52]. 

The spectra of each component were obtained from the fluorescence emission 
spectra (Figure 6). The [Eu(TODGA),]** (n = 1 — 3) complexes and Eu(III) solvent 
species were also discovered. No complexes with TODGA values greater than or equal 
to 4 were found in the fluorescence emission spectra, proving that the ratio of Eu(III) to 


TODGA equals 1:3. 
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Figure 5 Normalized Eu(III) emission spectra as a result of the transitions of 5Do—7F1 


and 5Do—7F2 for the complexation of Eu(III) with increasing TODGA in acetonitrile 


Eu** 
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Figure 6 Single-component spectra of the [Eu(TODGA):]** complexes (n = 1 — 3) and 
Eu(II) solvent species 

The slope value was disputable, although many studies showed that the Ln): 
DGA structures are 1:3 complexes. Antonio et al. elucidated the inner coordination 
sphere of Eu(III) and TODGA using X-ray absorption spectroscopy (XAS)[53]. One 
ether O atom and two carbonyl O atoms of a TODGA molecule make up each ligand’s 
tridentate to Eu(III). The coordination number is nine, with twelve distant carbon 
neighbors, six closest ones from carbonyl carbon atoms, and six slightly farther ones 
from the ether carbon atoms. Turanov et al.[54] found that the calculated slope value of 
Eu(iII) extracted by TODGA in n-decane containing 0.002 mol/L dinonylnaphtalene 
sulfonic acid (HDNNS) was 3.28 + 0.04. This was estimated by the plot of logDeu — 
log[TODGA] and the existence of EuLn(NOsz)3 species by the slope of logDeu — 
log[HNOs], which was equal to 2.86 + 0.08. Pathak et al.[55] investigated the influence 
of solvent type, acidity of extraction, and concentration of TODGA on luminescence 


lifetime and found that various ratios of TODGA to Eu ions resulted in the formation of 


Eu(TODGA)3** species when Eu(III) was complexed with TODGA. According to the 
findings of Sasaki et al.[56], extracted Eu(I1) complexes require three or four TODGA 
molecules and three nitrate ions to maintain stability in non-polar diluents by slope 
analysis. Zhu et al.[57] performed a slope analysis of TODGA to extract some Lns(III) 
and found that the slope was 3.9 for Eu(III), suggesting that the M(TODGA)4(NO3)3 
species finally formed. However, in this work, a 1:3 complexation ratio for TODGA 
with Eu(III) was determined. 

The theoretical Dg, values at different concentrations of TODGA can also be 
obtained using Equation (8). Equation (9) is obtained by linear fitting, as shown in 
Figure 4. 

LogDgy = 3.10 Log|TODGA] (org) + 6.75 (9) 

Equation (10) can be obtained from Figure 1 for irradiated HNO3-TODGA/nDD. It 
is an exponential radiolysis equation for the concentration of TODGA and the absorbed 
dose. 


CTODGA = -3 
r =e 2.20x107 °R (10) 
CTODGA 


We can obtain Equation (11) for Dg, and the absorbed dose by substituting 
Equation (10) into Equation (9). This empirical equation enables us to construct a 
mathematical relationship between the D;,, and absorbed dose in this system. 


log Dy, = —2.965 x 10-3R + 4.583 (11) 


@ Experimental Deu 
—— Theoretical D,,(Equation 11) 


logDe, = -2.965 x 10°°R + 4.583 
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Figure 7 Experimental log D;,, and theoretical log Dg, after irradiation for HNO3- 
TODGA/nDD 


The theoretical D;,,with absorbed doses ranging from 50 to 500 kGy was obtained 
using Equation (11). As shown in Figure 7, the experimental Dg, was higher than the 
theoretical D,;,, after irradiation, particularly at a high absorbed dose above 300 kGy 
for the HNO3-TODGA/nDD system. The experimental Dg, is close to the theoretical 
Dey, below 300 kGy, indicating that the decrease in Dgu is mainly related to the loss of 
TODGA concentration, with a slight influence from radiolytic products. 

Because some radiolytic products remained in coordination structures similar to 
those of TODGA, such as Pi and P2, which were identified by UPLC-QTOF, the 
complexation of radiolytic products on Eu(IIJ) was investigated further using DFT 


calculations. 


Theoretical calculations of complexation of TODGA and radiolytic 
products with Eu(IID) 


According to the slope and fluorescence titration analyses, the Eu(II): TODGA 
stoichiometry was 1:3. In addition, Kimberlin et al.[58] reported that several radiolytic 
products with the TODGA skeleton are involved in heteroleptic complexes with 
TODGA with a stoichiometry of 1:3 from the ESI-MS spectra. The time-resolved laser 
fluorescence results showed no change in the peak shape of the fluorescence spectra 
before and after irradiation (Figure S11). However, the fluorescence lifetime gradually 
decreased with increasing absorbed dose (Figure S12), indicating that y radiation did not 
change the M: L ratio of 1:3 but caused a slight partial change in the original fluorescent 
species. 

DFT calculations were employed to investigate the coordination abilities of 
TODGA and its radiolytic products with TODGA skeletons (Pı and P2). Figure 8 shows 
the optimized structures of the complexes formed by Eu(III) and the ligands. 
Theoretical calculations were performed for the [Eu(TODGA)2(P,)]** (n = 1 — 2) 
species and [Eu(TODGA)3]** complexes. It is realistic that one TODGA molecule in the 
initial complex is replaced by one radiolytic product to form mixed Ln-radiolytic 
product-TODGA complexes during y irradiation, which was already found in the ESI- 
MS spectra by Kimberlin et al.[58]. In addition, we investigated the optimized 
structures and coordination abilities of [Eu(Pn)3]** (n = 1 — 2), which were formed using 
three radiolytic products to represent a more extreme coordination environment of 


Eu(IID) after y irradiation. 
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Figure 8 Optimized structures of the ligands (TODGA, Pı and P2) and complexes 
formed by Eu(III); the ligands (light pink, blue, green, red, and white spheres represent 
Eu(iI), N, C, O, and H, respectively). 


The changes in the Gibbs free energy, entropy, and enthalpy for the complexes 
formed by the ligands (TODGA, Pı and P2) and Eu(III) are shown in Table 3. It can be 
seen that the formation of [Eu(TODGA)3]** complex has more negative AG in the gas 
phase (-2527.3 kJ/mol) and in the nDD phase (-1115.0 kJ/mol) than that of the 
[Eu(TODGA)2(Pn)|** and [Eu(Pn)3]** (n = 1 — 2) complexes in the two phases. 


Moreover, the coordination abilities of the radiolytic products of TODGA decreased, as 


reflected by the changes in the Eu—O bond length (Table 4). The average Eu — O bond 
length of the complexes changed from 2.4911 A ({Eu(TODGA)3]**) to 2.5000 A 
({Eu(P2)3]**), indicating a decline in the coordination ability of the complex. These 
results indicate that TODGA shows better coordination ability than the radiolytic 
products of TODGA; the complexation order is TODGA > Pi > P2. However, the 
formation of the [Eu(TODGA)2(Px)]** and [Eu(P»)3]** complexes (n = 1 — 2) had a high 
negative AG in the two phases. Theoretical calculations also indicated that the radiolytic 
products of TODGA, such as Pı and Pz, still maintain a good coordination ability with 
Eu atoms, proving that the radiolytic products retain partial complexation for Eu(III). 
Consequently, the experimental Dg, of TODGA was higher than the theoretical Dgu 
from the loss of TODGA with the increase in radiolytic products at a high absorbed 
dose. Furthermore, if the contents of Pi and P2 are determined exactly, by combining the 
TODGA contents, we can predict the Dz, of the irradiated TODGA/nDD system more 
precisely. This study will be conducted in future works. In our previous experiments, 
we found that the radiolytic products Pi and P2 were difficult to synthesize, and the 
radiolytic product Pi was not sufficiently stable. Therefore, we used DFT calculations to 
evaluate the coordination abilities of these radiolytic products, thereby justifying the 
experimental result that the experimental Dg, was higher than the theoretical Dg, due 


to the high dose TODGA losses in this study. 


Table 3 Changes in Gibbs free energy, entropy, and enthalpy (298.15 K, kJ/mol) for the 
complexes formed by ligands (TODGA, Pı and P2) and Eu(III) were acquired separately 
in (a) gas phase, as well as (b) nDD at the level of B3LYP/6-31G(d)/RECP. 


(a) Complexation in the gas phase AH, AGg TAS, 


Eu** + 3TODGA —> [Eu(TODGA)3]** -2693.7 -2527.3 -166.4 
Eu? + 2TODGA+ P; > [Eu(TODGA),(P)}** -2681.7 -2506.1 -175.6 
Eu** + 2TODGA+ P, —> [Eu(TODGA),(P2)}** -2657.8 -2487.8 -170.0 
Eu** + 3P; —> [Eu(P:)3]** -2674.4 -2462.8 -211.6 
Eu** + 3P2 —> [Eu(P2)3]** -2568.0 -2380.8 -187.2 
(b) Complexation in the nDD phase AH soi AG oi TASsgo1 
Eu** + 3TODGA —> [Eu(TODGA)3]** -1281.3 -1115.0 -166.3 
Eu** + 2TODGA+ P; — [Eu(TODGA))(P1)]}** -1276.4 -1100.9 -175.5 
Eu** + 2TODGA+ P; —> [Eu(TODGA),(P2)}** -1265.4 -1095.4 -170.0 
Eu** + 3P; > [Eu(P;)3]}** -1282.1 -1070.5 -211.6 
Eu** + 3P2 —> [Eu(P2)3]** -1226.4 -1039.2 -187.2 


Table 4 Changes in Eu — O bond length of complexes formed by Eu(IID) and ligands 
(unit: A) 


Structure Eu-O Eu-O Eu-O Eu-O- Eu-O Eu-O Eu-O- Eu-O- Eu-O Average 
[Eu(TODGA)3]3* T ta ao ie 2.436 2.439 oo ~ oF 2.4911 
[Eu(TODGA):(P1)] 2.421 2.422 2.426 2.430 2.436 2.453 2.602 2.611 2.627 2.4924 
a 9 2 6 6 3 2 4 1 0 : 
[Eu(TODGA):(P2)] 2.386 2.414 2.414 2.421 2.423 2.541 2.581 2.595 2.687 2.4964 
3+ 8 2 7 5 7 5 3 7 9 ` 
2.418 2.421 2.427 2.437 2.438 2.440 2.613 2.614 2.621 
3+ 
[Eu(P1)3] 7 7 0 1 1 8 8 6 3 2.4926 
2.362 2.382 2.393 2.479 2479 2.514 2.604 2.627 2.655 
3+ 
[Eu(P2)3] 8 1 1 8 9 7 3 1 9 2.5000 


In conclusion, by combining the exponential equation for the radiolysis kinetics of 
TODGA in nDD and the equation for slope analysis of the concentration of TODGA 
and Deu, we can derive a mathematical equation that includes the absorbed dose and Deu, 
which could be used to predict the extraction performance at low absorbed doses. This 


method can be used to study other extractant systems. This work provides in-depth 


insights into the effect of y irradiation on the extraction behavior of the TODGA/nDD 


system and provides useful information on the use of nuclear fuel reprocessing. 


4 Conclusions 


In summary, y-radiolysis of TODGA/nDD and HNO3-TODGA/nDD_ was 
investigated and compared with y-radiolysis of neat TODGA by analyzing the loss of 
TODGA and liquid radiolytic products. When the absorbed dose was increased, the 
TODGA concentration decreased exponentially, and the dose constant of TODGA in 
nDD was higher than that of neat TODGA because of the “sensitization effect” of nDD. 
However, the radiolysis of TODGA in nDD was only slightly affected by the pre- 
equilibration with HNOs:. Seven radiolytic products of TODGA were identified and 
semi-quantified using UPLC-QTOF-MS, and four radiolysis routes for TODGA were 
proposed. The Dgu of the irradiated systems decreased after irradiation, but the HNO3- 
TODGA/nDD system irradiated with an absorbed dose of 500 kGy maintained a high 
Dry of 2.2x10°. For the irradiated HNO3-TODGA/nDD, we obtained an empirical 
equation between D,,, and the absorbed dose by combining the radiolysis kinetics of 
TODGA and the complexation equation of TODGA with Eu(III), which fits well with 
the experimental results with absorbed doses below 300 kGy. DFT calculations 
demonstrated that the radiolytic products of TODGA with similar coordination 
structures possessed good coordination abilities with Eu atoms, which led to the 
experimental Dg, being higher than the theoretical Dg, based on the ligand content at 


high absorbed doses. 
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